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1. Introduction 

This paper pertains to a device for generating useful elec- 
trical power by means of an electrohydrodynamic (EHD) process. 

The device utilizes the flow of a fluid, normally a gas, such as 
air or water vapor, in which are entrained a very large number of 
very fine and well distributed solid or liquid particles, for 
example, water droplets as in an aerosol spray. The particles are 
electrically charged, either positively or negatively, as may be 
convenient in a particular case. The gas is caused to flow through 
a nozzle-like channel by the imposition of a suitable pressure drop. 
Let station 0 be the inlet to the nozzle. The charged particles 
are introduced into the stream of gas by a suitable injector at 
some upstream location, call it station i, and are removed from the 
stream by a suitable collector grid at some downstream location, 
call it station 2. An important intermediate location, station 1, 
is at the throat of the nozzle. The electrical charges which are 
transported downstream along with the particles on which they re- 
side constitute an electrical current. These charges move through 
an electrical field which exerts forces upon them in a direction 
and sense opposed to the general fluid motion. By proper design of 
the injector, the size of the water droplets or particles may be so 
regulated that these droplets have low mobility, that is, they move 
approximately with the surrounding gas with negligible relative slip 
produced by the opposing electrical forces. 

The gas stream does work on the charged particles in moving 
them against the resistance of the electrical forces. In this 
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process the gas stream undergoes a corresponding decrease in 
total enthalpy. The work done upon the charged particles creates 
a difference in electrical potential between stations i and 2. 

These stations are connected by an external circuit which includes 
a useful electrical load. Thus the enthalpy drop of the gas is 
ultimately converted into a useful electrical power output from 
the external circuit. 

The above scheme accomplishes the primary conversion of ther- 
modynamic energy directly into electrical power without the use 
of any major rotating or reciprocating components such as large 
turbines or electrical generators. Nevertheless, there still 
exists a requirement for a small pump to recirculate the condensate 
and a fan to recirculate the carrier fluid. A condenser and a 
boiler are also normally a part of the complete system, but these 
various auxiliaries are not considered here in any detail as our 
present concern is primarily with the nozzle in which the basic 
thermo-electric power conversion occurs. 

The above concepts are well known and can be found in the 
technical literature. See References 1, 2 and 3, for example. 

They are reviewed in this introduction merely to provide a proper 
background for our result which is quite specific and which is a 
consequence of the following analysis. 

The performance that can be achieved by an EHD device is 
limited by, among other things, the maximum electrical field 
strength that can be sustained at the most critical point in the 
field without inducing electrical breakdown, through a spark 
discharge, of the carrier medium. The electrical field that is 
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present in most EHD devices has its highest value at the point 
of injection. See References 1-4, for example. In this regard 
there is experimental evidence that breakdown field strength in 
any region varies linearly with the quantity Rp where R is 
the gas constant and where p is the fluid density in the region 
considered. This experimentally established limit depends of 
course on the type of gas/aerosol combination that is involved. 

But the presence of a maximum field at the beginning of the 
EHD conversion channel limits performance severely because the 
rest of the channel must operate below maximum capacity or, in 
other words, the electric pressure is at its allowable limit in 
only a small fraction of the conversion channel. 

The above considerations suggest that an optimum design 
would be one in which the local field strength is everywhere uni- 
formly close to the critical limit. Such a design would achieve 
greater electric work output per unit mass of fluid than would 
any other. The essence of this paper is the recognition that a 
uniform maximum loading through the channel is the optimum load- 
ing and that this condition can indeed be achieved at least to 
an acceptable approximation. A secondary aim of this paper is* 
to derive the characteristic geometrical features of the optimum 
design and certain associated performance parameters and limits. 

The geometrical, electrical, and thermodynamic features of 
such a power generator are governed by various physical laws of 
which one of the most significant in the present context is 
Poisson's equation as it applies to an electric field. In its 
full generality, Poisson's equation is three dimensional in 
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nature; by restricting attention to configurations having polar 
symmetry we can simplify this to a two-dimensional form. More- 
over, by further restricting the application to configurations 
whose largest radial dimensions are small in comparisons with 
their axial length, we may finally simplify Poisson's equation to 
a one-dimensional approximation. Despite certain limitations in 
accuracy, this one-dimensional version can be more enlightening 
than the more elaborate two-dimensional analysis because it shows 
basic trends so much more simply and clearly. It is this one- 
dimensional approximation on which the present paper is based. 

The analysis must deal not only with Poisson's equation for 
electrical fields, but also with the laws of fluid flow. It is 
consistent with the foregoing commentary to write these laws also 
in their common one-dimensional forms. At this level of ideal- 
ization it is also appropriate to treat the flow through the 
channel as isentropic. We also consider the fluid as a perfect 
gas for which y , the ratio of specific heats, is constant. 

One other special assumption is involved. This is based 
on the fact that the drop in total temperature through the 
channel is normally very small in comparison with the absolute 
total temperature T q at the nozzle inlet, station 0. Hence it 
is. - permissible for certain analytical purposes to neglect this 
temperature drop and to approximate the flow through the channel 
as an isentropic flow of constant total temperature. Calculations 
to be presented later show that the actual change of total tem- 
perature amounts to less than one percent for a typical case. 
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Future progress in EHD power generation might increase 
electric power output by an order of magnitude over levels that 
seem currently feasible. In that case the drop in total tempe- 
rature through an EHD nozzle could become significant. Fortunate- 
ly, the present analysis can be revised and generalized to account 
rigorously for such variation in total temperature when and if 
necessary. Under present circumstances, however, it is not 
warranted to complicate the analysis to include this refinement 
since at low power outputs it has no appreciable effect on the 
final calculated results. 

2. Analysis 

It is convenient for the purposes of the present analysis to 
introduce a parameter a which is used to distinguish between 
the two distinct cases of positively and negatively charged par- 
ticles. Specifically, we set a = + 1 for the case of positively 

charged particles, and we set o = - 1 for the case of negatively 

charged particles. It is then appropriate to denote the electrical 
charge per unit mass by means of the product aq where q has 
units of coulombs/kg and is always positive by definition while 
a is a dimensionless factor as defined previously. 

Using the above notation we may write the following two 
expressions for the electric current i and the mass flow rate 
m through the channel, namely, 

i = cqpAv = constant along channel (2.1) 

ih = pAv = constant along channel (2.2) 
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where 



i = electric current, amps 

ih = mass flow rate through channel, kg/sec 
q = electric charge per unit mass, positive by 



definition, coulomb/kg 

p = fluid density, kg/m^ (variable) 

A = cross-sectional area of channel at any 

2 

given streamwise station, m (variable) 
v = mean axial velocity at any given streamwise 
station, m/sec (variable) 

Note that Eq. (2.1) is based on the previously stated 
assumption that size of the fluid particles and the charge on 
each particle can be so regulated that the mobility of the par- 
ticle with respect to the surrounding medium is essentially 
negligible . 

By dividing Eq . (2.1) by Eq. (2.2) we also find that 

t = crq = constant along channel (2.3) 

m v 

The one-dimensional version of Poisson's equation can now 
be written in the form 




where z = axial coordinate, m 

$ = electrical potential, volts 
e = = electrical permittivity of free space 

= 8.854 x 10 -2 farad/m 



(2.4) 
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The local electrical field strength at an arbitray axial 
station z is defined as usual 



E = - 




(2.5) 



where E = field strength, volts/m 

As will be seen here, an optimum EHD generator has an elec- 
tric field which does not change sign from inlet to outlet. 

Under these circumstances the electric field is always negative 
for positive space charged and positive for negative space charge. 
Hence we may write 



E = - c|E| (2.6) 

It has been established by experiment that over a broad (but 
not unlimited) range, the field strength at breakdown is well 
represented by the simple linear law (see Reference 2) . 




= C o + C B R P 



(2.7) 



where C and C_ are characteristic constants of the medium, 
o B 

R is the gas constant in units of Joule/kg°K . Experimental 
measurements show that constants C q and Cg happen to have the 
same numerical values for both air and steam. However, we defer 
to later pages of this paper any reference to actual experimental 
values. Of course, the present analysis is restricted to that 
specific range of conditions for which Eq. (2.7) is in fact a 
valid approximation. 

If we now impose the constraint that the local field strength 
is everywhere just equal to its critical value at impending break- 
down, we readily find from the last three equations that 
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HI - °< c o + s r p > 



( 2 . 8 ) 



Now differentiating this gives 




= a c„ R fi 

B dz 



(2.9) 



Upon equating the right sides of Eqs. (2.4) and (2.9) and 
simplifying, we obtain the important relation 



I ($£.) = _ q 

p I dz I e C B R 



( 2 . 10 ) 



The reciprocal of the quantity on the right side of Eq. (2.10) 
(momentarily disregarding the negative sign) now identifies a 
significant characteristic length; let us denote it by symbol X . 

£ C R 

X = — — (meters) (2.11) 

q 



It is very instructive to rewrite Eq. (2.11) in the follow- 
ing alternative form 



Aq “ A . = £ C R = constant 

mm max B 



( 2 . 12 ) 



Notice that the quantity z C R is a characteristic pro- 
perty of the medium. Once the medium is chosen the designer has 
no further control of the value of this constant. Eq. (2.12) 
now tells us that in order to maintain the electrical loading at 
incipient breakdown, the product Xq must remain constant. Thus 
any increase in q must be accompanied by a corresponding decrease 
in X , or vice versa. 
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Experience up to now indicates, however, that there exists 

some practical upper limit ( 3 inax to the value of q that can 

actually be achieved for any specific type of design, and hence 

there exists some corresponding lower limit X . on the accom- 

panying longitudinal characteristic length. Of course, a value 

of q less than q max can always be employed if necessary in 

which case the corresponding value of X will be greater than 

X . . 

mm 

The situation is shown schematically in Fig. 2.1 in the form 
of a log-log plot. The solid line represents the locus of points 
all of which produce incipient breakdown. Incidentally, it will 
also be shown later that all of these points correspond to the 
same theoretical power output from the generator. Note that the 
line continues on indefinitely toward decreasing values of q 
and increasing values of X . At the right end, however, it ter- 
minates at the point corresponding to q and X . . 

max mm 

Some typical values of these various characteristic constants 
based on published data are summarized in Table 2.1 for reference 
purposes . 

Upon substituting X into Eq. (2.10) , we can easily inte- 
grate the simple expression that is thereby obtained. The result 
is 




(2.13) 
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X (log scale) 




q 



( log scale ) 



Fig. 2.1 Design Points and Design Limits 

for an Optimum Slender EHD Nozzle 
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Table 2.1 Typical EHD Properties 



Source of Data 


kef. 1 


Ref. 2 


Ref. 3 


Medium 




Water 
droplets 
in steam 


Water 
droplets 
in air 


Water 
droplets 
in air 


Quantity 


Units 


C 

o 


volts/m 


8.63 x 10 5 


8.63 x 10 5 


8.63 x 10 5 


C B 


m 2 °K/C 


9.49 x 10 3 


9.49 x 10 3 


9.49 x 10 3 


£ C 0 R = Aq 


C m/kg 


3.9 x 10 -5 


2.4 x 10 -5 


2.4 x 10" 5 


g max 


C/kg 


2.5 x 10 -3 


3 x 10 -2 


3 x 10~ 2 


A . 
min 


m 


1.5 x 10" 2 


8 x 10 -4 


8 x 10~ 4 



This important result is the key to all remaining details of 
the analysis. 

Eq. (2.13) reveals that the density continues to drop mono- 
tonically with increasing distance downstream. Thus the flow cor- 
responds to the known isentropic flow through a converging-diverging 
nozzle for which the walls are so designed that the density drops 
exponentially as specified by Eq. (2.13) . Let subscript 1 denote 
the throat of the nozzle. Then the following relations are known 
to apply 




(2.14) 
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/ \ 




1 + Y M 2 1 




1 = 


2 


KJ 


1 


Y + 1 

2 J 



(Y + 1) 

W'-'TT 



1 

M 



(2.15) 



where y 
M 
r 

r l 
It is 



= ratio of specific heats 
= Mach number at station z 
= nozzle-Tadius at station 
= nozzle radius at throat 
readily deduced from Eqs. 




(Y - 1) 



In 



1 + 



Y 



z 



(2.13) and 2.14) that 




(2.16) 



Eqs. (2.15) and (2.16) now constitute a pair of parametric 

Z 

equations for ( — ) as a function of ( T ) with Mach number M 

r l A 

as parameter. These equations define the shape of the optimum 
nozzle having uniform loading just below breakdown all along its 

IT Z 

length. A tabulation and a curve of ( — ) versus (y) are 

r l A 

given in the next section, using X = 1.3 . This value of X is 
representative for steam. See Table 4.1 and Fig. 4.1. 

At the throat of the nozzle M = M^ = 1 . Then Eq. (2.16) 

gives 



(Y - 1) 



£n 



Y + 1 



For y = 1.3 this yields 



(2.17) 




(2.18) 
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Let z^ denote the axial station which locates the injector 



that introduces unipolar electrical charges into the flow. We 
do not stipulate in advance any particular fixed location for the 
injector because the nature of the charged particles may require 
one location or another depending on certain supersaturation con- 
ditions that might be needed for the formation of the aerosol 
droplets (see Reference 1) . Moreover, the converging-diverging 
nozzle does not have to be of the shape specified in this analysis 
upstream of the injector location. 

Of great importance is the theoretical gross power output of 
such an optimum device. To determine this we must first find the 
current from Eq. (2.1). It is advantageous to evaluate the terms 
on the right side of Eq. (2.1) at the nozzle throat, station 1. 
Thus , 



Now letting a Q and a^ denote sonic velocity at stations 
0 and 1, respectively, we have 



Upon combining the above two relations we finally obtain the 
current in the form 



i = aqp^v 



(2.19) 




( 2 . 20 ) 




( 2 . 21 ) 
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Note, however, that from the previous analysis of the isen- 



tropic nozzle flow, it follows that 




( 2 . 22 ) 



Moreover, it is useful to eliminate q from Eq. (2.21) by 
means of Eq. (2.12) and to rearrange the result in dimensionless 
form. In this way we finally obtain the dimensionless current in 
the form 



(Y + 1) 

2(y - 1) 

(2.23) 

The next step is to find the change in electric potential 
between stations i and 2. This is accomplished by integrating 
Eq. (2.8) between these limits. We again rearrange the result 
in dimensionless form. The result is 



/ Ai \ — cr 


2 


C B R ^o a o A l ) 


1 

* — 1 

+ 

1 




(2.24) 



Finally, the theoretical power output follows from tije rela- 
tion 



P e = i(<0 2 - 40 



(2.25) 



where = theoretical gross electrical power, watts 
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